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METHODS AND APPARATUS FOR ANALYSIS OF 
ANGIOGRAPHIC AND OTHER CYCLICAL 
IMAGES 

BACKGROUND OF THE INVENTION 

[0001] This invention relates generally to methods and apparatus for 
diagnostic medical imaging and more particularly to methods and apparatus for 
automated angiographic analysis. However, the invention is also more generally 
useful for analysis of other cyclical images. 

[0002] A coronary angiogram (or synonymously, an arteriogram) is 
an x-ray of the coronary arteries. Such images help the physician observe blocked 
arteries for the diagnosis of coronary artery disease. An angiogram is often part of a 
cardiac catheterization procedure including other tests and measurements. A patient 
undergoing catheterization also undergoes physiological monitoring. Typically, an 
electrocardiogram (ECG or, synonymously, EKG) is produced from electrical leads 
attached to a patient and analyzed to detect irregularities indicative of heart damage or 
disease. Data from the ECG is stored in a computer database. 

[0003] To perform an angiogram, a physician injects a contrast 
medium into the coronary arteries using catheter inserted through a blood vessel (such 
as a femoral artery) up to the heart and then takes an X-ray of the heart. The 
angiogram assists the physician in determining the location and severity of coronary 
artery disease. This information is then used to determine what, if any, treatment is 
medically appropriate. 

[0004] Cardiac angiography is performed in a cardiac catheterization 
laboratory utilizing cardiovascular X-ray imaging equipment. Left ventricular 
analysis, which is the contour detection of the diastolic and systolic left ventricle, can 
be carried out by capturing images in which a contrast medium has been injected into 
the left ventricle. 
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[0005] Images produced by X-ray imaging equipment can be stored 
on photographic film or electronically as digital images. Typically, digital images are 
stored in DICOM (Digital Imaging and Communications in Medicine format) and 
stored and viewed at a workstation. 

[0006] Cardiac angiography images are usually reviewed as moving 
pictures. For left ventricular analysis, the perimeter of the dye filled left ventricular 
chamber is manually or semi-automatically delineated in each portion of the cardiac 
cycle. Various known analysis methods require a user to manually select a pair of 
images, more particularly, a diastole and a systole image, draw or mark the contours 
of the dye filled ventricle. Automatic (i.e., computerized) edge detection is possible, 
but is often difficult due to lack of edge definition or obscuration by musculature 
inside the ventricle. More often than not, automated edge detection systems require 
manual adjustment for successfiil edge detection for left ventricular analysis. 

BRIEF DESCRIPTION OF THE INVENTION 

[0007] Some configurations of the present invention therefore 
provide a method for analyzing stored firames of time-stamped imaging data and time- 
stamped physiologic data obtained during a study of a patient during a plurality of 
cardiac cycles. The imaging data includes firames representative of a succession of 
images of the patient's heart. These configurations use a computer to analyze the 
stored physiologic data to determine times of a selected phase in the plurality of 
cardiac cycles, select, in the plurality of cardiac cycles, stored images corresponding to 
the determined times of the selected phase, and align and average images 
corresponding to the selected phase across the plurality of cardiac cycles to produce a 
single representative frame of image data. 

[0008] Some configurations of the present invention provide a 
cardiology analysis apparatus for analyzing stored frames of time-stamped imaging 
data and time-stamped physiologic data obtained during a study of a patient during a 
plurality of cardiac cycles. The imaging data includes fi:ames representative of a 
succession of images of the patient's heart. In the apparatus, a computer is configured 
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to analyze the stored physiologic data to determine times of a selected phase in the 
plurality of cardiac cycles, select, in the plurality of cardiac cycles, stored images 
corresponding to the determined times of the selected phase, and align and average 
images corresponding to the selected phase across the plurality of cardiac cycles to 
produce a single representative frame of image data. 

[0009] Still other configurations of the present invention provide a 
method for analyzing stored frames of time-stamped cyclical image data obtained 
during a study of an object during a plurality of cycles. The method uses a computer 
to analyze the stored time-stamped cychcal image data to determine time of a selected 
phase of the cycle, select, in the pliu-ality of cycles of time-stamped cyclical image 
data, stored images corresponding to the determined times of the selected phase, and 
align and average images corresponding to the selected phase across the plurality of 
cycles to produce a single representative frame of image data. 

[0010] Still other configurations of the present invention provide an 
analysis apparatus for analyzing stored frames of time-stamped cyclical imaging data 
and time-stamped waveform data obtained during a study of an object during a 
plurality of cycles. The imaging data includes frames representative of a succession of 
images of the object. The apparatus includes a computer that is configured to analyze 
the stored waveform data to determine times of a selected phase in the plurality of 
cycles, select, in the plurality of cycles, stored images corresponding to the determined 
times of the selected phase, and align and average images corresponding to the 
selected phase across the plurality of cycles to produce a single representative frame 
of image data. 

[0011] It will be appreciated that configurations of the present 
invention directed to cardiac image analysis provide fiilly automatic left ventricle 
analysis using multiple images of systole and diastole. Using these multiple images, a 
single representative image of each phase proves more accurate and repeatable edge 
detection and dependent measures, such as ejection fraction than is known in current 
methods and apparatus. Corresponding advantages can be obtained in configurations 
of the present invention directed to other imaging applications, including non-medical 
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applications in which various forms of cyclical data and images of an object are 
recorded. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] Figure 1 is a block diagram depicting a network for 
facilitating storage and review of data acquired in a cardiac catheterization laboratory. 

[0013] Figure 2 is a block diagram depicting a cardiac catheterization 
laboratory networked to an asynchronous DICOM server with NTP time 
synchronization server. 

[0014] Figure 3 is a flowchart indicating data acquisition by the 
hemodynamic monitoring system depicted in Figure 2. 

[0015] Figure 4 is a flowchart indicating data acquisition by the x-ray 
imaging system depicted in Figure 2. 

[0016] Figure 5 is a flowchart indicating data processing by the 
workstation of Figure 1. 

DETAILED DESCRIPTION OF THE INVENTION 

[0017] A technical effect of the systems and processes described 
herein include at least one of producing a single representative frame of image data 
from a plvirality of stored images, providing left ventricle analysis using multiple 
images of systole and diastol, and facilitating measurement of ejection fraction, 

[0018] In some configurations of the present invention and referring 
to Figure 1, a local area network (LAN) 18 facilitates communication between a 
systems housed in a cardiac catheterization laboratory 12, a physician review or 
overview workstation 14 and a DICOM server 16. For example, angiographic x-ray 
images acquired by imaging equipment in laboratory 12 and formatted as DICOM 
objects are transmitted by LAN 18 to a database 19 accessed by DICOM server 16. 
Workstation 14 can be used by a physician to display and perform quantitative 
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analysis on retrieved images. Certain analysis, including left ventricular analysis, 
require the selection of images corresponding to particular cardiac events (e.g., systole 
and diastole). Workstation 14, in many configurations, is a computer that is part or all 
of a cardiology analysis apparatus. The cardiology analysis apparatus may further 
comprise an operator interface 13, such as a mouse and/or a keyboard, and a display 
monitor 15 and/or printer 17. 

[0019] Referring now to Figure 2, cardiac catheterization is 
perfomied on a patient 26 in a catherization laboratory 12 while patient 26 is resting 
on an examination table 28. After patient 26 is anesthetized, a catheter is placed in the 
heart of patient 26 by inserting a catheter into blood vessels in the neck, arm, or groin, 
and advancing the catheter through the blood vessels to the heart. The catheter is 
maneuvered in the heart, and a contrast agent is injected when the catheter is in 
position while a rapid sequence of X-ray images are recorded by X-ray imager 30. X- 
ray imager 30 is controlled by a technician located at X-ray control station 52. X-ray 
images can be viewed on live video monitors 32 and 46 located in catheterization 
laboratory 12 and control room 20, respectively, hi an angiography procedure that 
provides a roadmap of blood vessels, these live X-ray images show blood flowing 
through the heart chambers or through blood vessels connected to the heart. A 
contrast-enhanced roadmap image that aids in catheter guidance is displayed on 
roadmap monitors 34 and 48 located in catheterization laboratory 12 and control room 
20, respectively. 

[0020] Each frame of imaging data acquired by X-ray imager 30 is 
time-stamped. The time is determined by reference to a first clock 31 located in or 
associated with X-ray imager 30. Each frame is later encapsulated in a DICOM 
object, with the time stamp associated with that frame inserted in a predetermined 
header field. The DICOM object from X-ray imager 30 can then be uploaded to 
DICOM server 38 via DICOM image network 24 for storage in a DICOM image 
database 19 (see Figure 1). DICOM server 38 is programmed with database 
management software configured to manage the DICOM image database. The 
DICOM image database can reside on DICOM server 38 itself or elsewhere on 
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another server, computer, or storage system in communication with DICOM server 
38. 

[0021] Cardiac catheterization laboratory 12 also houses a 
hemodynamic monitoring system 40, for example, a MacLab 7000 system 
manufactured and sold by GE Medical Systems Mformation Technologies, Mc, 
Milwaukee, WI. This system is operated by a technician who is seated at a 
hemodynamic control station 50 in control room 20. During catheterization 
procedures, the patient undergoes physiological monitoring. Hemodynamic 
monitoring system 40 records physiologic data, including ECG waveforms acquired 
fiom electrodes placed on the patient. Data acquired by hemodynamic monitoring 
system 40 is displayed on hemodynamic monitors 42 and 44 located in catheterization 
laboratory 12 and control room 20, respectively. Hemodynamic monitoring system 40 
is also programmed to perform ECG analysis, producing signals (i.e., "pointers") that 
indicate features of the ECG waveform, such as exact diastolic and systoUc heart 
beats. Hemodynamic monitoring system 40 can also automatically acquire and 
display 12-lead ECGs, continuous ST segment analysis, respiration rate, 
thermodilution cardiac output, Sp02 level, noninvasive blood pressure, and up to four 
invasive pressures. 

[0022] Hemodynamic monitoring system 40 also maintains a catheter 
laboratory event log ("cath lab event log") in which the various procedures being 
performed are recorded. The hemodynamic technician controls all information in the 
cath lab event log. This technician has his or her hands free while watching and 
listening to the entire catheter laboratory procedure behind leaded glass. The 
technician enters the projection being shot and the artery selected for live video 
imaging into the cath lab event log before the injection of dye. 

[0023] The pointers indicative of features of the ECG waveform are 
time-stamped by hemodynamic monitoring system 40 utilizing a second local clock 
41 located in or associated with hemodynamic monitoring system 40. The time- 
stamped ECG waveforms are later encapsulated in DICOM objects with the 
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associated time stamps in predetemiined header fields. These DICOM objects are 
uploaded to DICOM server 38 via DICOM waveform network 22. 

[0024] Quantitative imaging analysis is performed on acquired X-ray 
images using software modules resident on X-ray control station 52, on hemodynamic 
control station 50 or on any physician review workstation capable of communicating 
with DICOM server 38, such as station 14 represented in Figure 1. Such analysis 
requires the selection of images acquired substantially in synchronization with certain 
cardiac events of interest. This selection could be accomplished by selecting firames 
of imaging data having a time stamp closest in time to time-stamped pointers 
corresponding to particular cardiac events. However, first clock 31 and second clock 
41 must be synchronized to avoid time offset differences (or, synonymously, 
"offsets"). Such offsets may be significant and may change in magnitude over time. 

[0025] To synchronize first clock 31 and second clock 41, a Network 
Time Protocol (NTP) time synchronization server software module 51 resides in 
DICOM server 38. In addition, data acquisition systems 30 and 40 and clocks 31 and 
41, respectively, are also each provided with respective NTP time synchronization that 
synchronize clocks 31 and 41 of data acquisition systems 30 and 40. Each NTP time 
synchronization client software module initiates a time request exchange with the 
NTP time synchronization server software module 51. As a result of this exchange, 
the client is able to calculate the link delay, its local offset, and adjust its local clock to 
match a reference clock at another computer (e.g., DICOM server 38) where server 
module 51 resides. NTP time synchronization modules for clocks 31 and 41 also 
update their respective clocks (for example, at regular one minute intervals) to ensure 
that they do not drift significantly fi"om synchronization. 

[0026] Accordingly, firames of imaging data acquired by X-ray 
imager 30 can be synchronized with features in the physiological data acquired by 
hemodynamic monitoring system 40 that represent predetermined cardiac events of 
interest. For example, fi-ames of imaging data representing a succession of images of 
a patient's heart during a study are acquired using imager 30. Imager 30 is 
programmed to time stamp each acquired firsune with the a time stamp utilizing clock 
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31. Imager 30 is also programmed to determine respective offsets of clock 31 relative 
to a reference clock such using, for example, NTP client time synchronization module 
for clock 31 and NTP time synchronization server software module 51. Acquired 
frames and time stamps and local offsets associated with the acquired frames are 
uploaded to DICOM server 38 via DICOM image network 24. This data is 
transmitted in the format of DICOM objects, each DICOM object encapsulating a 
respective frame of imaging data and having a header comprising respective fields for 
a time stamp and an offset associated with that frame. 

[0027] During the same study, physiologic data is acquired from the 
patient using hemodynamic monitoring system 40. Hemodynamic monitoring system 

40 is programmed to time stamp acquired physiologic data corresponding to one or 
more predetermined cardiac events. The time stamps are generated utilizing clock 41. 
Hemodynamic monitoring system 40 is also programmed to determine offsets of clock 

41 relative to the reference clock used by NTP time synchronization server 51 for each 
time stamp associated with the physiologic data. The acquired physiologic data with 
associated time stamps and local offsets are uploaded to DICOM server 38 via 
DICOM waveform network 22. This data is transmitted in the format of DICOM 
objects, each of which encapsulates physiologic data and which includes a header 
comprising fields for a time stamp and an offset associated with that physiologic data. 

[0028] DICOM server 38 executes a database management software 
module that instructs DICOM server 38 to store imaging and physiologic DICOM 
objects in a database (not shown in Figure 2). This database management software 
module also instructs DICOM server 38 to respond to requests for files from remote 
workstations, hi particular, DICOM server 38 can be requested to send a set of 
DICOM files to a physician review station for quantitative imaging analysis. A 
quantitative analysis software module can be resident, for example, on x-ray control 
station 52, hemodynamic control station 50 or physician review station 14 (see 
Figures 1 and 2). The quantitative analysis software module interprets time stamps 
and local offsets in DICOM object headers and merges images with physiologic data 
(e.g., ECG waveforms) for automatic image selection. For example, in one 
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configuration, when a physician replays images, all diastole frames and matching 
systole fi-ames are displayed. The displayed firames are selectable for left ventricular 
analysis. Individual fi-ames for stenosis evaluation are identified as to their 
contractility in the DICOM header. The quantitative analysis software module allows 
a physician to select an individual frame for quantitative coronary analysis (QCA) 
with an indicator to indicate whether or not it is at the peak when the physician selects 
a QCA tool. The physician may then frame forward or backward until both the 
required position and contractility is attained. 

[0029] In some configurations, review station 14 is a cardiology 
analysis system having an operator interface, a display monitor and a computer. The 
computer is programmed to compensate for unsynchronized local clocks of imager 30 
and hemodynamic monitoring system 40. Compensation is accompUshed using the 
associated local offsets when time stamps of frames of imaging data are compared to 
time stamps of physiological data corresponding to cardiac events of mterest. More 
particularly, to obtain the true times of data acquisition relative to the reference clock 
used by NTP time synchronization server 51, the computer is programmed to adjust 
the time stamps by the amount of offset before the time stamps are compared. M this 
manner, the quantitative imaging analysis software module is able to select a frame 
that was acquired substantially in synchronization with the occurrence of a 
predetermined cardiac event recorded in the physiological data for the study in issue. 
The selected frame of imaging data is then utilized for quantitative coronary analysis. 
For example, two frames corresponding to the diastole and systole of a particular 
cardiac cycle are selected and then analyzed for left ventricular analysis. 

[0030] In yet another configuration, local clocks 31 and 41 at imager 
30 and hemodynamic monitoring system 40, respectively, are synchronized with the 
reference clock used by NTP time synchronization server 51 in real-time using NTP 
time synchronization. By "real time," it is meant that clocks 31 and 41 are adjusted so 
that their time stamps represent adjusted, i.e., synchronized, tune. In this 
configuration, the header of each DICOM object has a field containing the adjusted 
time stamp, so that a local offset field is not required. 



-9- 



133865 



Express Mail No. EL 977935812 US 



[0031] Configurations of the present invention are not limited to 
using NTP. Other distributed network time protocols that are used in various 
configurations include, for example. Simple Network Time Protocol ("SNTP"), 
Experimental Network Time Protocol ("XNTP"), Secure Network Time Protocol 
("STIME"), and Real-Time Transport Protocol ("RTP"). 

[0032] Figure 3 is a flowchart representing a physiological 
monitoring process that is performed in catheterization laboratory 12 in some 
configurations of the present invention. Leads connected to a patient undergoing the 
procedure pick up electrical pulses to and firom the patient's heart at 54. The 
waveforms of the electrical pulses are sampled and displayed on a monitor as 
waveforms at 56 and stored in a computer memory at 60. Hemodynamic monitoring 
system 40 is programmed to perform ECG analysis of waveforms at 58, results of 
which are also stored in the computer memory. As part of that analysis, hemodynamic 
monitoring system 40 is programmed to time stamp pointers to various characteristic 
features in the acquired waveforms. Optionally, quantitative analysis (QVA volume 
measurement 66 or QCA size measurement 68) can be performed at hemodynamic 
control station 50. QCA and QVA software modules may be included as part of 
hemodynamic monitoring system 40, or one or both of these software modules may 
reside on a PC connected to hemodynamic monitoring system 40. While the 
procedure is under way, a hemodynamic technician enters identifying codes and 
activates time stamps for various procedures being performed at 62. These procedures 
are maintained in a procedure log 64 called the "cath lab event log". 

[0033] hi some configurations and referring to Figure 4, various steps 
of the imaging process are performed in catheterization laboratory 12. X-ray images 
are acquired in some configurations at a fi^e rate of 8 to 60 fi-ames per second at 72. 
Acquired images are displayed as live video at 74 on live video monitor 32 (see 
Figure 2). A contrast-enhanced roadmap image is displayed at 76 on roadmap 
monitor 34 (see Figure 2). The roadmap image assists a physician guiding the catheter 
through the patient's blood vessels and can also be used to perform quantitative image 
analysis (e.g., QCA or QVA) at 78. hnages are time stamped and stored in an image 
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archive 80, along with, in some configurations, results fi-om quantitative image 
analysis. Image archive 80 contains the database managed by DICOM server 38 (see 
Figure 2). DICOM server 38 also dovmloads DICOM data from image archive 80 to a 
requesting workstation at 82, for example via a TCP/IP network, hnages can then 
undergo quantitative image analysis as described above. 

[0034] In some configurations of the present invention, an 
angiographic LV (Left Ventricular) analysis is performed utilizing multiple dye-filled 
images at diastole and systole. An edge detection algorithm using multiple images 
facilitates total image edge detection of dye-filled left ventricles, providing enhanced 
accuracy and repeatability. In these configurations, images of the dye filled ventncle 
at the same phase of different cardiac cycles provide additional information than 
would be obtained by analyzing a single image at that phase. 

[0035] For example, Figure 5 is a fiowchart indicating data 
processing by workstation 14 of Figure 1, in which the technical effect is achieved by 
analyzing, from a DICOM database of a patient study, electrocardiograph (ECG or 
synonymously, EKG) waveforms with acquisition time stamps 102 to find systolic 
and diastolic times at 104. Also from the DICOM database of the patient study, image 
frames of the left ventricle 106 are selected at 108. The selected images are diastolic 
and systolic images based on ECG images and acquisition timestamps. More 
particularly, images from two or more cardiac cycles are selected at 108, so that there 
are two or more times and frames for both phases. 

[0036] For the pluraUty of selected systolic images, the images are 
aUgned with one another at 110 utilizing a aUgmnent software module and averaged at 
112 utilizing an averaging software module to produce a single representative systolic 
image (An "image," as that term is used in Figure 5 and in comiection therewith, 
more accurately refers to a digital representation of an image.) Portions of the image 
that do not represent left ventricular anatomy are filtered out or removed from the 
single image at 1 14 utiUzing a filtering software module. The edge of the venfricular 
wall is then detected at 1 16 utilizing an edge detection software module, and a systolic 
volume is determined at 118 from the detected edge utilizing a volume-determimng 
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software module. A corresponding sequence of operations 120, 122, 124, 126, and 
128 is also performed on diastolic images to determine a diastolic volume. From both 
the determined diastolic volume and the determined systolic volume, some 
configurations of the present invention determine an ejection fraction utilizing an 
ejection fraction determining software module. 

[0037] Thus, some configurations of the present invention provide a 
method for analyzing stored frames of time-stamped imaging data 106 and time- 
stamped physiologic data 102 obtained during a study of a patient during a plurality of 
cardiac cycles. The physiologic data may be ECG waveform data, pressure 
measurements, or other data. A computer, such as workstation 14, is used to analyze 
104 the stored physiologic data 102 to determine time of a selected phase of the 
plurality of cardiac cycles. The stored data in some configurations is X-ray imaging 
data, which may have been enhanced by injection of a contrast agent into the 
cardiovascular system of patient 28. Computer 14 is also used to select 108, in the 
plurality of cardiac cycles, stored images 106 corresponding to the determined times 
of the selected phase and to align 110 and average 112 images corresponding to the 
selected phase across the pluraUty of cardiac cycles to produce a single representative 
frame of image data. 

[0038] In some configurations, computer or workstation 14 is further 
used to filter 114 non-left ventricle anatomy from the single representative frame of 
image data. Also in some configurations, computer or workstation 14 is used to detect 
1 16 edges of a left ventricle wall in the filtered, single representative frame of image 
data, and is fiirther used in some configurations to determine 118 cardiac volume in 
the filtered, single representative fiame of image data. 

[0039] M some configurations, the selected cardiac phase is a systoUc 
phase, and the determined cardiac volume is a systolic volume, hi various other 
configurations, the selected cardiac phase is a diastoUc phase and the determined 
cardiac volume is a diastohc volume. In still other configurations, both a selected 
systoUc phase and a selected diastolic phase are processed to determine both a systolic 
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volume and a diastolic volume, and computer 14 is used to determine an ejection 
fraction utilizing the determined systolic volume and the determined diastolic volume. 

[0040] Some configurations of the present invention provide a 
computer or workstation 14 configured to analyze the stored physiologic data to 
determine times of a selected phase in the plurality of cardiac cycles. The computer is 
further configured to perform all or some subset of the various processes and 
procedures described above. 

[0041] Configurations of the present invention are useful in 
applications other than cardiac imaging. For example, some configurations utiHze 
time-stamped cyclical waveform data having recognizable phases along with a 
plurality of related cycles of time-stamped cyclical image data collected during a study 
of an object. In these configurations, the time-stamped cycUcal waveform data is used 
in much the same way as the EKG or ECG waveform data, except that systolic and/or 
diastolic phases are replaced by other phases appropriate to the intended application. 
Also, images for the intended application other than those of a patient's heart are used. 
In these applications, images corresponding to a selected phase in different cycles are 
averaged to produce a single image representative of the object being imaged at the 
selected cycle. Filtering may also be performed that corresponds to non-left 
ventricular anatomy filtering in a cardiac application, and edge detection suitable for 
the application may also be performed. If the object being imaged is a pumping 
apparatus, volumes at two phases can be determined analogously to the determination 
of the systolic and diastolic volumes, and an ejection fi-action (or corresponding 
measure) may be determined. 

[0042] It will thus be appreciated that various configurations of the 
present invention provide fully automatic left ventricle analysis using multiple images 
of systole and diastole. Using these multiple images, a single representative image of 
each phase proves more accurate and repeatable edge detection and dependent 
measures, such as ejection fraction than is known in current methods and apparatus. 
Corresponding advantages can be obtained in configurations of the present invention 
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in other applications, including non-medical applications in which various forms of 
cyclical data and images of an object are recorded. 

[0043] While tiie invention has been described in terms of various 
specific embodiments, those skilled in tiie art will recognize that the invention can be 
practiced witii modification within tiie spirit and scope of ttie claims. 
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